Complex viscosity behavior and cluster formation in attractive colloidal systems 
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The increase of the viscosity, which is observed in attractive colloidal systems by varying the 
temperature or the volume fraction, can be related to the formation of structures due to particle 
aggregation. In particular we have studied the non trivial dependence of the viscosity from the 
temperature and the volume fraction in the copolymer-micellar system L64. The comparison of 
the experimental data with the results of numerical simulations in a simple model for gelation 
phenomena suggests that this intriguing behavior can be explained in terms of cluster formation 
and that this picture can be quite generally extended to other attractive colloidal systems. 

PACS numbers: 83.80.Uv, 82.70.Dd, 83.60.Hc 
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Soft materials (polymers, proteins, colloids, etc.) rep- 
resent a research field of large interest in science and tech- 
nology and are a common domain for physics, chemistry 
and biology. The effective interaction amongst particles 
in these systems is in many cases the result of the compe- 
tition between repulsive and attractive contributions and 
may give rise to aggregation processes with the formation 
of a variety of meta-stable or stable mesoscopic structures 
This is the case of gelation phenomena in at- 
tractive colloidal systems at low volume fractions 
Here slow dynamics and eventually structural arrest are 
associated to the formation of long-living non-compact 
structures, with a change in the viscoelastic properties 
as the one typically observed in gelation phenomena. At 
present, this is a point of central interest in the study of 
attractive colloidal systems as the effects of aggregation 
processes at low volume fractions on the dynamics are 
still poorly understood 0- |T^|. 

A first interesting clue comes from the better known gel 
formation in polymer systems: here the phenomenology 
of gelation is typically interpreted in terms of a percola- 
tion transition |13|- 115| . that is the formation of a per- 
manent spanning structure, accompanied by a structural 
arrest. This suggests that gelation phenomena could be 
understood in terms of formation of clusters, and even- 
tually of a spanning structure, also in attractive colloidal 
systems. Here, on the other hand, at high volume frac- 
tions and low temperatures, structural arrest phenomena 
are observed with typical glassy features, that have been 
interpreted within the framework of the mode coupling 
theory 0,0,0]. In this case the glassy behavior can be 
strongly influenced by the short range attraction, but it 
cannot be simply described in terms of formation of clus- 
ters or of a percolating structure. Therefore the role of 
structure formation in the slow dynamics and the struc- 
tural arrest of attractive colloidal systems is not com- 



pletely clear. 

In this work we discuss the viscosity behavior in the 
copolymer micellar system L64, where the short range 
attraction is due to an effective intermicellar interaction. 
We observe a transient gelation phenomenon followed by 
a structural arrest. For a fixed value of the temperature, 
as function of the volume fraction the viscosity increases, 
reaches a plateau and then increases again at higher vol- 
ume fraction. In this process the rheological properties 
exhibit a crossover from gelation to a structural arrest 
typical of the glass transition. The comparison of the ex- 
perimental findings with the numerical study of a simple 
model suggests that these two phenomena, which appear 
to be quite distinct, stem from two different mechanism: 
One related to the presence of long living clusters, the 
other due to crowding of particles. The model also sug- 
gests that the interrupted gel phenomena can be more 
generally related to colloidal gelation In other col- 

loidal systems, due to different features of the effective 
attractive interaction, the gelation phenomenon occurs 
at much lower temperature and volume fraction and it 
is more pronounced. Consequently, structural arrest and 
gelation almost coincide and occur close to the percola- 
tion transition. At higher temperatures clusters can be 
less persistent and structural arrest, occurring at higher 
volume fractions, exhibits features of the so called at- 
tractive glasses. Increasing further the temperature the 
effect of the bonding interaction vanishes and one finds 
features of the hard sphere glass transition 0, 0, 0] . 

We have studied the viscosity behavior, as a function 
of the temperature T and of the micellar volume frac- 
tion <!>, in the copolymer-micellar system L64 [l9|-[23|. 
This is an aqueous solution of a non-ionic three-block 
copolymer made of polyethylene and polypropylene ox- 
ides, [{PEO) 13 - {PPO) 3Q - (PEO) 13 ] (L64 for short). 
These polymer molecules in water are T dependent sur- 
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factants that spontaneously form hydrated monodisperse 
spherical micelles in a wide T — $ range. At low tem- 
peratures, both PEO and PPO are hydrophilic, so that 
the L64 chains readily dissolve in water, and the polymer 
molecules exist as unimers. However, increasing T, PPO 
tends to become less hydrophilic than PEO creating an 
unbalance in the hydrophilicity between the central block 
and the end block of the polymer molecule that conse- 
quently acquires surfactant properties and self-assemble 
to form micelles. The inter-particle attraction is mainly 
due to the fact that, on further increasing T, water be- 
comes progressively a poor solvent to both PEO and 
PPO chains. The phase diagram of the L6A/D2O system 
is shown in Fig. ^ Due to the intermicellar attraction, 
the system has an inverted binodal curve with a lower 
consolute critical point, a T — $ dependent percolation 
line where the viscosity exhibits a steep increase, and 
a glassy line as predicted by the mode coupling theory 
|l9l l2f]j. In particular, by increasing the volume fraction 
the percolation line precedes the glassy line. 
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FIG. 1: The phase diagram of the L64/D2O system (adapted 
from Ref. |20||): the inverted binodal line (circles), the per- 
colation line (empty squares) and the glassy line (circles and 
triangles on the right of the figure). The percolation and the 
glassy lines were measured by both quasi elastic light scatter- 
ing and small angle neutron scattering. 

Viscosity measurements, as a function of T at differ- 
ent <&, were performed with a strain-controlled rheomctcr 
using a double wall Couette geometry, at the fixed fre- 
quency lu = 1 sec^" 1 . To ensure a linear response, for 
different T and <&, we have tested our system as a func- 
tion of low applied strains, 7. After this check a low 
strain deformation, 7 = 0.05, was maintained for all the 
experiments. A crucial point is that the Peclet number 



Pe — 7£ 2 /2Z?o, which characterizes the amount of distor- 
tion of structures with linear dimension £ (Dq is the par- 
ticle short-time diffusion coefficient at low concentration 
where hydrodynamic effects dominate) must be small. 

In Ref. H3 it was found that £ - 500ol and D Q of the 
order of 2 • 10e~ n m 2 /s close to the percolation thresh- 
old, ensuring Pe < 1. These distortions are responsible 
for shear thinning in the system structures, which is ex- 
pected to set on at Pe ps 1. In other words, the condition 
Pe < 1 constitutes an important experimental constraint 
in the study of the elasticity of systems characterized by 
clustering processes like polymer solutions and attractive 
colloids. In particular, this is true for viscosity measure- 
ments in phenomena dominated by long-range "critical" 
correlations, like those which originate the diverging in- 
cipient cluster typical of the percolation threshold in the 
sol-gel transition. Just in these processes, shear thinning 
can originate a "levelling off" of the measured viscosities. 

In Fig. [2 the measured viscosity, 77, of the system is 
plotted as a function of the temperature, T, at differ- 
ent volume fractions, $. At each <I>, 77 is character- 
ized by a growth from the values typical of the solvent 
(77 ~ 10 _2 Pa) followed by a bending at high tempera- 
tures. At high <j> a second growth is observed in 77 by 
further increasing the temperature: The first growth of 
the viscosity corresponds to the percolation transition 
(PT) 0, and the second growth to the structural arrest 
(SA) of the attractive glass transition |22| . 
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FIG. 2: The L64/ D2O viscosity, 77, as a function of the tem- 
perature, T, for different values of the volume fraction $. In 
figure the percolation transition (PT) and the structural ar- 
rest (SA) are shown. 

In Fig. |2| 77 is plotted as a function of the volume frac- 
tion, $, for different values of the temperature, T. At low 
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FIG. 3: The L64/ D2O viscosity, 77, as a function of the vol- 
ume fraction, $, for different values of the temperatures, T, 
in the range 10 — 46° C. The continuous line is a power law, 
($ c - $)~ fe , with k = 1.35. 

volume fraction a first strong growth of the viscosity is 
found close to the percolation line. However, a bending 
in the viscosity is found, as the volume fraction is in- 
creased above the percolation threshold, and finally, by 
further increasing the volume fraction, a second growth 
corresponding to the glass transition 22] is observed. A 
similar effect in the viscosity behavior has also been ob- 
served in another micellar system in Ref. [24|. 

The strong growth of the viscosity close to the perco- 
lation line is reminiscent of the viscosity divergence ob- 
served in polymer gelation, and might be due to the pres- 
ence of long living clusters in the system. The intermicel- 
lar interaction has an entropic origin. In fact, the micelle 
can be schematized as a hard sphere with polymers at- 
tached on the surface. When two micelles come close, 
the polymers get entangled (long- living bonds). Due to 
the entanglement, the two micelles have to overcome an 
entropic barrier to get separated. Therefore the bond 
lifetime t& will be related to the height of such entropic 
barrier Sb, i.e. Tf, ~ e Sb (the entropy is in units of Boltz- 
mann's constant). The other quantity of interest is the 
characteristic time t u for two close micelles to get entan- 
gled (bonded): r u ~ e " where S u is the entropic barrier 
needed to be overcome to go from the un-bonded state 
to a bonded state. Therefore the probability pb for two 
close micelles to be bonded is given by 

P b e Sb _|_ e S„ 

Note that such bond probability can be quite smaller 



than 1, depending on the value of S u . This could explain 
why in micellar systems percolation can occur at quite 
high volume fraction, as compared to other colloidal sys- 
tems 

To study the effect of aggregation and structural ar- 
rest in such micellar system, we propose a model which 
contains these two characteristic times for forming and 
breaking bonds. Tuning the parameters Tf, and t u we 
can study the gelation process with infinite bond lifetime 
(polymer gelation) and finite bond lifetime (micelles). A 
detailed discussion on the model and on the numerical 
simulations can be found in Ref. |25| . 

We consider a solution of tetrafunctional monomers 
with excluded volume interactions and performed Monte 
Carlo simulations on the cubic lattice. At t = the vol- 
ume fraction $ is fixed, and bonds between monomers 
are randomly quenched. The four possible bonds per 
monomers, randomly selected, are formed along lattice 
directions between monomers that are nearest neighbors 
and next nearest neighbors. In the case of polymer 
gelation, once formed the bonds are permanent. The 
monomer diffuse according to the bond-fluctuation dy- 
namics (BFD)[26(. We let the monomer diffuse to reach 
the stationary state and then study the system for dif- 
ferent values of the concentration. The system has a 
percolation transition at $ c = 0.376 ± 0.003. In the case 
of colloidal gelation, we introduce a finite bond lifetime 
Tf, . We start with the same configurations of the previous 
case, with a fixed $ where the bonds have been randomly 
quenched. During the monomer diffusion with BFD at 
each time step we attempt to break each bond with a fre- 
quency 1/Vf,. Between monomers separated by a distance 
less than l bonds are then formed with a frequency 1/t„. 

From the time autocorrelation function, f^{t) , of equi- 
librium density fluctuations with q= (f > -j j f )i we calcu- 
late the relaxation time r as the time such that fq(r) ~ 
0.1. In Fig. 0] the relaxation time, t, is plotted as func- 
tion of the volume fraction, <I>, for the permanent bonds 
and for the finite lifetime bonds case at different values 
of Tf,. In the figure one finds the data for the permanent 
bond case on the left, and from left to right the data for 
finite bond lifetime, for decreasing values of T5. 

In the case of permanent bonds, t($) displays a power 
law divergence, ($ c — $) _fc , at the percolation thresh- 
old, with k = 1.35. For finite bond lifetime Tf, the re- 
laxation time instead increases following the permanent 
bond case, up to some value $* and then deviates from 
it. The longer the bond lifetime the higher 4>* is. For 
higher $ the increase of the relaxation time corresponds 
to the onset of the glassy regime in the relaxation behav- 
ior discussed in Ref. [25j. This qualitative behavior has 
also been confirmed via molecular dynamics simulations 
in Ref. [2~3]. This truncated critical behavior of the re- 
laxation time qualitatively reproduces the experimental 
data shown in Fig. [3J In order to stress the analogy, we 
have plotted a power law behavior for the viscosity with 
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FIG. 4: The average relaxation time, r, as a function of the 
volume fraction, $. From left to right: the data for the per- 
manent bond case diverge at the percolation threshold with 
a power law (continuous line), see text. The other data refer 
to finite r b = 3000, 1000, 400, 100 MC step /particle decreas- 
ing from left to right (the dashed lines are guides to the eye). 
(Adapted from Ref. 



the same value of k = 1.35 (continuous curve in Fig. 

This comparison suggests the following interpretation 
for the experimental data. We consider that clusters of 
different sizes will be present in the system at any time 
t. In the permanent bond cluster of radius R 

diffuses amongst the others with a characteristic relax- 
ation time, t(R); At the percolation threshold the con- 
nectedness length critically grows in the system and so 
does the overall relaxation time. In the case of a finite 
bond lifetime, Tb, there will be a cluster size R* so that 
Tb ~ t(R*). That is, clusters of size R < R* can be con- 
sidered as permanent, whereas clusters of size R > R* 
will not contribute with their full size to the enhance- 
ment of the relaxation time in the system. Therefore, a 
long but finite bond lifetime produces the enhancement of 
the viscosity due to the formation of persistent clusters. 
On the other hand, it introduces an effective cluster size 
distribution with a cut-off which keeps the macroscopic 
viscosity finite in the system [2^ . 

In conclusions, based on a minimal model for colloidal 
system we have analyzed the experimental viscosity data 
in copolymer-micellar systems. For a fixed temperature 
the viscosity exhibits a first steep growth as the volume 
fraction increases. This effect is due to coalescing clus- 
ters of monomers linked by long living bonds. Thus the 
viscosity increases as in a typical gelation phenomenon, 
but due to the finite lifetime of the bonds this increase is 



interrupted and followed by a saturation. As the volume 
fraction increases further, the viscosity exhibits a second 
steep growth until structural arrest is reached. The re- 
sulting glass transition can be related to the formation of 
a spanning cluster made of localized particles connected 
by bonds. It is the presence of the bonds which makes 
the difference between the attractive glass and a repulsive 
glass. 
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